Momentum spectra, anisotropic flow, and ideal fluids 
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If the matter produced in ultrarelativistic heavy-ion collisions reaches thermal equilibrium, its 
subsequent evolution follows the laws of ideal fluid dynamics. We show that general predictions can 
be made on this basis alone, irrespective of the details of the hydrodynamical model. We derive 
several scaling rules for momentum spectra and anisotropic flow (in particular the elliptic flow, v%, 
and the hexadecupole flow, V4) of identified particles. Comparison with existing data is briefly 
discussed, and qualitative predictions are made for LHC. 
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An ultrarelativistic Au-Au collision at RHIC produces 
a dense system of interacting particles and fields, which 
then expands into the vacuum. Ideal-fluid models have 
been succesful in describing this expansion they are 
able, to a certain extent, to reproduce the magnitude and 
the transverse momentum (p t ) dependence of elliptic flow 
of identified particles, together with their p t spectra, for 
momenta pt ^ 2GeV/c. 

In this Letter, we derive general properties of momen- 
tum spectra of identified particles emitted by an ideal 
fluid, which do not depend on the specific model used. 
We shall introduce an important distinction between slow 
particles, whose velocity equals the fluid velocity at some 
point, and fast particles, whose velocity exceeds the max- 
imum fluid velocity. We discuss in detail the implications 
of ideal-fluid behaviour for both slow and fast particles. 

The expansion of a plasma into a vacuum goes through 
successive steps. If the particle mean free path inside the 
plasma is small enough, the plasma thermalizes, leading 
to a collision-dominated, isentropic expansion, which is 
described by ideal (i.e., inviscid) fluid dynamics. On the 
contrary, the late stage of the expansion is collision-free. 
In between, a transition regime takes place. This transi- 
tion regime must in principle be modelled by transport 
theory, which encompasses both isentropic and collision- 
free limits HQ. 

In the context of heavy- ion collisions, however, the 
transition regime is most often modelled by a simple 
Ansatz, the sudden freeze-out approximation, which is 
a sharp transition between the two extremes: one first 
defines a space-time hypersurface E along the history of 
the ideal fluid, on which the transition is expected to 
take place. At each point of E, free-streaming particles 
are emitted according to thermal distributions in the rest 
frame of the fluid. Integrating over E, one obtains for a 
given particle the following momentum spectrum jj]: 
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where u^(x) is the fluid 4- velocity at point 1 on E, C 



is a normalization constant, and we have neglected the 
effects of quantum statistics (in practice, the latter may 
only be significant for pions at lowpt). For simplicity, we 
also assume that the fluid temperature T is everywhere 
the same on E, but the results derived in this paper do 
not rely on this assumption. The possibility has also 
been raised that particles of different types |5| or with 
different transverse momenta pt [f| have different freeze- 
out temperatures. 

Let us briefly discuss the validity of the Cooper-Frye 
Ansatz, Eq. lfT)l. This approximation is expected to be 
poor when applied to observables which are sensitive to 
the detailed physics at freeze-out, such as HBT radii [tJ. 
The reason is that Eq. assumes an isotropic momen- 
tum distribution in the rest frame of the fluid, while 
non-relativistic studies have shown that freeze-out pre- 
cisely occurs when the relative difference between par- 
allel and transverse components of the kinetic tempera- 
ture becomes large ||- On the other hand, the simple 
freeze-out Ansatz may be a reasonable one for comput- 
ing single-particle spectra, provided that collective ex- 
pansion dominates over random, thermal motion. This 
is known as the hypersonic approximation in the context 
of non-relativistic gas dynamics 0,0- In the language of 
heavy-ion collisions, it can be rephrased as follows: con- 
sistency of the ideal-fluid picture requires that the freeze- 
out temperature be much smaller than the inverse slope 
parameters obtained by exponential fits to transverse- 
mass spectra 0. The best-fit value of the freeze-out 
temperature at RHIC is T - 100 MeV HJ03, while the 
inverse slope parameter for pions is T e g = 210 MeV [l3j . 
The value of T is large enough to expect significant devia- 
tions from ideal- fluid behaviour, i.e., viscous effects [l4j . 
It is however interesting to study the small-T limit in 
view of upcoming heavy-ion experiments at LHC, and 
also to have a better grasp on viscous effects, which are 
easily seen as deviations from this limit. 

We therefore investigate systematically the properties 
of momentum spectra in the limit of small T. The general 
idea is that the integral over E in Eq. can then be 
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FIG. 1: (Color online) Schematic representation of the dis- 
tribution of the fluid transverse 4-velocity. The :r-axis is the 
impact-parameter direction. The short (resp. long) arrow in- 
dicates the four-velocity of a slow (resp. fast) particle. The 
most probable values of the fluid four- velocity for this particle 
are marked as darker (red) areas. 



performed by means of a saddle-point integration. 1 In 
physical terms, it means that the dominant contribution 
comes from the points where the energy of the particle 
in the fluid frame, p^Uy., is minimum. For a given p^, 
p^Uy is a function of the space components of the fluid 4- 
velocity, u (the fourth one being related to them through 
u° = \J 1 + u 2 ) , which themselves depend on the point x 
on S. Since p^u^ is the energy of the particle in the fluid 
rest frame, its absolute minimum is the particle mass m. 
This minimum is reached when the particle is at rest with 
respect to the fluid, i.e., when its velocity p/p° equals the 
fluid velocity u/u° (or, equivalently, if u — p/m). 

This absolute minimum, however, occurs only if there 
exists a point on £ where this value of the fluid velocity 
is reached. This leads us to a qualitative discussion of the 
values taken by u at freeze-out. The longitudinal fluid ve- 
locity is expected to span almost the whole range from — 1 
to 1 in ultrarelativistic collisions, and the simple Bjorken 
picture u z /u° = z/t 0] shows, at least qualitatively, 
how it is related to space-time coordinates. The radial 
fluid velocity spans a more limited range. The reason is 
that transverse collective flow is not initially present in 
the system but builds up progressively. For a given fluid 
rapidity Uf = h ln((u° + u z ) / (u° — u z )), the transverse 4- 

velocity, ^Ju 2 , + u 2 , extends up to some maximum value 
u maxi which may depend on yf, and on the azimuthal 
angle <j)f for non-central collisions. As shown schemati- 
cally in Fig. [2 u ma x ( <f>f) is largest at 4> f = 0, along 
the direction of impact parameter. This is due to larger 
pressure gradients in this direction [l7| , which explain the 
large in-plane elliptic flow observed at RHIC [lq. Typi- 



cal values of u max are of order 1 at RHIC (u max = sinh p 
in the notations of with the best-fit value po ~ 0.9). 

From now on, we make a distinction between "slow" 
and "fast" particles as follows: a particle of mass m, with 
rapidity y and transverse momentum p t , is defined as 
slow if pt/m < u m ax(2/, 0) for all <j> (that is, actually, for 
(j) = 7r/2, where the minimum occurs). Conversely, a fast 
particle is defined by pt/m > u mavi (y , (p) for all <j> (that 
is, for = 0, where the velocity is maximal). Between 
both regimes, there is a small intermediate region, which 
will not be considered in this Letter. 

For a slow particle, there is a point on E such that the 
fluid velocity equals the particle velocity, and the mini- 
mum p^Uy — m is reached. (Our results are unchanged if 
there are two or more such points.) If T is small enough, 
the dominant contribution to the integral Q comes from 
the neighbourhood of this point (see Fig.^). The integral 
can then be evaluated approximately by expanding the 
exponent to second order around the minimum of p^u^. 
What remains is a Gaussian integral. For a given veloc- 
ity, p^ is proportional to the particle mass m. Hence, the 
width of the Gaussian varies with m like 1 / y/rh, and the 
integral over Y, in Eq. is mT^/ 2 times a function of the 
particle velocity. This means that the mass dependence 
is only a global factor for slow particles: 
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where / is the same for all particles. As a result, trans- 
verse momentum and rapidity spectra (integrated over 0) 
of identified slow particles coincide, up to a normalization 
factor, when they are plotted as a function oipt/m and 
y. The coefficients quantifying azimuthal anisotropics 
v n = (cos 710), which are independent of the total yield, 
should coincide for different identified slow particles at 
the same pt/m and y. This property holds for directed 
flow v\ as well as for elliptic flow vi and the higher har- 
monics. The behaviour is qualitatively correct for t>2 at 
RHIC, which rises more slowly with p t for heavier parti- 
cles nim. 

The condition under which the saddle-point approxi- 
mation is good for slow particles can be roughly stated 
as T < m for a relativistic fluid, for dimensional rea- 
sons: the larger m, the smaller the width of the Gaussian, 
and the better the approximation. A detailed calculation 
gives the condition 

T < m< ax (3) 



1 The same method was used earli er t o predict the be- 
haviour of longitudinal HBT radii Il5l . 



with v max = M max /w max and u max = \Jl + u max ; this 
condition amounts to assuming that collective motion 
dominates over thermal motion. At RHIC, we expect 
the approximation to be poor for pions (furthermore, the 
pion spectrum is contaminated at low pt by secondary de- 
cays, and may also be sensitive to Bose-Einstein statis- 
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tics), but it might be a reasonable one for kaons and 
heavier hadrons with p t < m/c. 

Let us now discuss fast particles. Roughly speaking, 
these are the particles that move faster than the fluid: 
the minimum value of p^u^ is larger than m. In order to 
locate this minimum, we denote by y (resp. yf) the par- 
ticle (resp. fluid) longitudinal rapidity, by ut\ the trans- 
verse component of u parallel to the particle transverse 
momentum p t , and by u± the transverse component of 
u orthogonal to p t . With these notations, 

p^Ufi = m t cosh(y - yj)Jl + u| + u\ - p t u\\ , (4) 



p|. Minimization with respect to 



where m t = \frn 
yf and uj_ gives yf — y and uj_ = 0, i.e., the fluid 
velocity is parallel to the particle velocity. The mini- 
mum is then attained when tiii is maximum, i.e, when 
Mil = Umax(y,<p), where <p is the azimuthal angle of the 
particle. In other words, fast particles come from regions 
on S where the parallel velocity is close to its maximum 
value (see Fig.^J. A saddle-point integration then gives 2 



dN 



dyd 2 p t y/pl 



: exp 



T 



(5) 

where the (y, <p) dependence is implicit. This result 
was already obtained long ago for massless particles in 
Ref. UK (see also 

The saddle-point approximation is valid for fast par- 
ticles if p t is large enough. A more precise criterion is 



m t u u 
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(6) 



together with the condition p t > mu max . At RHIC, 
the left-hand side (lhs) of Eq. © is smaller than the 
right-hand side (rhs) by at least a factor of 2 as soon as 
pt > 0.7 GeV/c for pions, pt > 1.2 GeV/c for kaons, and 
pt > 1.8 GeV/c for (anti)protons. On the other hand, 
ideal fluid dynamics is expected to break down if pt is 
too high, since high p t particles have been shown to be 
more sensitive to off-equilibrium (viscosity) effects [l4^. 
Deviations from fluid-like behaviour are best seen on el- 
liptic flow, for mesons above 1.5 GeV/c, and for baryons 
above 2.5 GeV/c. The window in which our approxi- 
mation works is likely to be narrow, which reflects the 
importance of viscous effects at RHIC. Better agreement 
should be reached at LHC. 

The pt spectra of identified particles are directly ob- 
tained from Eq. neglecting the <fi dependence of u m ax- 



We assume that the maximum value u ma x is reached at an inner 
point of E. If it occurs at the edge of the fluid, there is no square 
root in the pre-exponential factor. 



Radial flow results in flatter m t -spectra for heavier par- 
ticles. In addition, Eq. J5J implies a breakdown of rat- 
scaling: the slope of the spectrum decreases with increas- 
ing m t for pions, and increases for protons, in qualitative 
agreement with experimental findings [13J. 

For non-central collisions, we can also obtain the 
anisotropic flow coefficients. We expand u max ((f>) in 
Fourier series, and neglect odd harmonics: 

{<t>) =u max (l + 2y 2 cos(2</>) + 2V 4 cos(4^) + ---). (7) 

The parameter V2 is of the order of 4 % for semi- 
central Au-Au collisions at RHIC. It is related to the 
parameter P2 of blast wave parameterizations [ill |23| 
by V% = P2/(2w m ax)- The (p distribution is obtained by 
inserting Eq. into (JSJ. If T is small enough, the <j> 
dependence in Eq. J5J is dominated by the exponential. 
Expanding the latter to first order in V2, one obtains 
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m t v n 



c). 



(8) 



A similar equation was already obtained in Ref. j 2 ■ -J | in 
the framework of a simplified fluid model, and was shown 
to fit RHIC data rather well. In particular, Eq. JHJ shows 
that the "mass ordering" which follows from Eq. |J5J for 
slow particles persists at high p t in hydro: at a given p tl 
heavier particles have smaller i>2- 

We finally make predictions for the hexadecupole flow, 
V4. We expand the exponential of Eq. j£| and look for 
terms in cos(40). To leading order one obtains two terms: 



v iVPt) = 7^ (Pt - Wt^max) 



2T 2 
VlUmi 



(p t - m t v max ) , 



(9) 



For large enough p t , the first term dominates over the 
second, which gives the simple, universal relation 



Vi{Pt) 



V2(Pt 



(10) 



Let us derive the domain of validity of this approxima- 
tion. If u max ((f>) is a smooth function of 0, one gener- 
ally expects V4 to be of order (V2) 2 . The condition for 
Eq. JTj| is then 



(Pt 



m t v n 



c). 



(11) 



At RHIC, the lhs is smaller than the rhs by at least a 
factor of 2 for pions with p t > 0.8 GeV/c. For heavier 
particles, Eq. © supersedes Eq. ljTT)l . 

Our result, Eq. is in contradiction with the state- 
ment that V4 is a sensitive probe of initial conditions |2-l| : 
on the contrary, we find a universal result, which can be 
directly used as a probe of ideal fluid behaviour, not of 
initial conditions. The experimental value found by the 
STAR Collaboration [25| is a factor of 2 to 3 higher than 



4 




P t [GeV/c] 



FIG. 2: (color online) Numerical results from ideal hydro- 
dynamics with Bjorken longitudinal expansion [lrj| . and a 
black-body equation of state. Initial conditions mimic a mid- 
central (6=8 fm) Au-Au collision at RHIC. The initial den- 
sity has been fixed to reproduce (pt) = 450 MeV/c, as mea- 
sured for pions at RHIC jl3l . Top: Solid line: dN/dyd 2 p t ; 
dash-dotted line, fit using Eq. JSJ. Middle: V2(pt)- Bottom: 
V4(pt)/v2 (pt) 2 . 

our prediction. Deviations from ideal-fluid behaviour are 
generally expected to yield higher values of v± . 

Our results for fast particles, Eqs. JSJ, JSJ, (fTOft . are 
compared to results from a numerical 3-d hydrodynami- 
cal calculation in Fig. [2 The calculation has been pushed 
to very large times, so that the small-T limit applies. The 
value of V4/V2 does not go exactly to 0.5 at large p t but 
rather to 0.63. This is due to the fact that the initial ec- 
centricity is large for this value of the impact parameter, 
and Eq. I|10|) is obtained through a leading order expan- 
sion in the anisotropy. We have checked numerically that 
agreement is better for lower values of b, where the ec- 
centricity is smaller. 

Before we come to our conclusions, let us compare our 
approach with the popular blast-wave one. The blast- 
wave parameterization, in its simplest form, assumes a 
unique radial velocity for the fluid |27]]; this framework 
has recently been refined to take into account the az- 
imuthal dependence of the fluid velocity |2^| and of the 
freeze-out surface [2^] in non-central collisions, and even 
a distribution of fluid velocities A few parame- 

ters (typically four) are then fitted to experimental data. 
Some of the results we derived above were already ob- 
tained within the blast-wave approach, namely the mass- 
ordering of the V2{pt) of different types. However, our 
present framework is more general in the sense that we 
do not assume a given fluid- velocity profile, but also more 
specific in the sense that we assume that collective motion 
dominates over thermal (random) motion. In addition, 
blast-wave fits treat slow and fast particles on an equal 
footing, ignoring the distinction between both types of 
particles. Although fitting the whole spectrum with a 



single formula is admittedly more convenient, it misses 
an important feature of the underlying physics, since slow 
and fast particles originate from different regions of the 
expanding fluid. In particular, fits using our formulas for 
fast particles may yield values of T and u max which dif- 
fer from blast- wave fits. Finally, our formulas are signifi- 
cantly simpler than blast-wave parameterizations, which 
involve special functions. 

We have obtained the following results for momentum 
spectra and anisotropics in the framework of ideal-fluid 
models using a saddle-point approximation of the mo- 
mentum distribution: 

• At low pt, identified particles of different masses 
have the same momentum spectra and anisotropics 
(up to a normalization for the spectra) , when plot- 
ted as a function of velocity variables y and Pt/m. 
This defines "slow" particles. This scaling is due 
to the fact that slow particles move with the fluid: 
they come from the regions where the fluid veloc- 
ity equals their velocity. The scaling is expected 
to be poor for pions. It is expected to break down 
when pt/m exceeds w ma x, the maximum value of 
the transverse 4-velocity of the fluid. u max may in 
general depend on the rapidity y, and reflects the 
underlying equation of state of the expanding mat- 
ter. 

• Fast particles, defined by Eq. (|HJ), all originate 
from the region where the fluid is fastest along 
the direction of the particle velocity. As a result, 
their transverse momentum spectra and azimuthal 
anisotropics at a given rapidity are uniquely deter- 
mined by three parameters it ma x, T, and V2, and 
given by Eqs. JSJl, ®, 013) ■ Comparing the v 2 of 
different particles should directly give the precise 
value of itmax, while transverse momentum spectra 
yield T. 

These results can be used as signatures of hydrodynamic 
evolution in heavy-ion collisions, and also as consistency 
checks of numerical ideal-fluid calculations. Ideal-fluid 
evolution leads to different behaviours for slow and fast 
particles. Some of the results obtained for fast particles 
(in particular for elliptic flow) are already known from 
blast-wave approaches. We have shown that they are in 
fact more general. The scaling rules for slow particles, 
which are evidenced here for the first time, should be 
further tested on available RHIC data. We expect all 
our results to be in closer agreement with data at LHC 
than at RHIC. In particular, we predict that the value of 
the ratio va/{v-i) 2 should be lower at LHC than at RHIC. 
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